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The pumping conductance of a disordered two-dimensional Chern insulator scales with increasing 
size and fixed disorder strength to sharp plateau transitions at well-defined energies between ordinary 
and quantum Hall insulators. When the disorder strength is scaled to zero as system size increases, 
the “metallic” regime of fluctuating Chern numbers can extend over the whole band. A simple 
argument leads to a sort of weighted equipartition of Chern number over minibands in a finite 
system with periodic boundary conditions: even though there must be strong fluctuations between 
disorder realizations, the mean Chern number at a given energy is determined by the clean Berry 
curvature distribution expected from the intrinsic anomalous Hall effect formula for metals. This 
estimate is compared to numerical results using recently developed operator algebra methods, and 
indeed the dominant variation of average Chern number is explained by the intrinsic anomalous 
Hall effect. A mathematical appendix provides more precise definitions and a model for the full 
distribution of Chern numbers. 
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I. INTRODUCTION 


. The problem of how conduction electrons move when the host 
3 §etal is magnetic dates back at least to the measurement by 
flail in 1879 of the effect that bears his name. The term “anoma- 
*0us Hall effect” refers to the observation that in many magnetic 
' Materials an electric field generates a transverse current density, 


jy Try U/: 


( 1 . 1 ) 


even in zero magentic field. Current belief is that in most ma- 
r4«rials the primary effect on conduction electrons of the mag¬ 
netism, even in a ferromagnet, comes from the modification of 
■nrhe electron band structure via spin-orbit coupling rather than 
am magnetic fields. Recent years have seen a revival of inter- 
£t including a great deal of progress on both the microscopic 
understanding of the effect 2 and its strength in various impor¬ 
tant materials^' these accomplishments are summarized in a 
-cd ccnt comprehensive review by Nagaosa et a/P 

One reason why this effect is particularly interesting is that it 
, believed to receive an “intrinsic” contribution from the Berry 
;phas^ of the band structure, independent of the detailed na- 
. ^jjre of the disorder. The Berry phase is a geometric property 
^at probes how the Bloch states of the band electrons evolve as 
Sthe crystal momentum changes. As first discussed by Karplus 
< -£thd Luttingei®, the Berry phase appears in the semiclassical 
equations of motion of an electron wavepacket, and this was 
argued to explain the anomalous Hall effect. Even for non¬ 
interacting electrons, the connection between the Berry phase 
and observed transport phenomena is quite subtle because scat¬ 
tering processes give additional contributions to the transverse 
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conductivity a xy . The goal of this paper is to use recent ad¬ 
vances in the understanding of the Berry phase in insulators, 
specifically the integer quantum Hall effect 

Tie 2 

axv = ~h' n e ( L2 ) 

to give a different picture of the metallic anomalous Hall con¬ 
ductivity. 

For most of this paper we will focus on two spatial dimensions, 
where (in the absence of spin-orbit coupling) the ultimate fate of 
a zero-temperature system is to be an insulator, either ordinary 
(Try = 0) or quantum Hall (Eq. |I.2[ ) . Our goal will be to show 
that the metallic anomalous Hall effect can still be understood 
by interpreting it as a disorder average over quantized insula¬ 
tors. In each disorder realization, the insulator is characterized 
by a set of topological invariants determining n in (Eq. |I.2[ ) , and 
a xy is effectively quantized. The average over disorder gives a 
non-quantized a xy , and we give an heuristic argument for why 
this averaged conductance should be determined in the weak- 
disorder limit by the intrinsic anomalous Hall formula. This 
approach gives a numerically accessible limit where the dom¬ 
inant, and possibly only, contribution to cr xy is the intrinsic 
contribution. We then apply a recently developed operator al¬ 
gebra methocP to compute the disorder-averaged conductance 
numerically for several models with nontrivial structure in a xy 
as a function of electron density, which adds support to the 
theoretical picture. 

We now introduce some basic notions and then describe the 
results more precisely. The integer n is well known to be quan¬ 
tized in two-dimensional insulators: it is determined by the sum 
of the Chern numbers or TKNN integer^ of the occupied bands. 
The Chern number is a topological invariant that is expressed 
as an integral over a band of the Berry curvature: defining the 
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Berry connection or vector potential as 

A J = KI-*VK>, (1.3) 

where j is a band index and k is crystal momentum, the Berry 
curvature is 


P = V fc x A- 7 ', 


(1.4) 


and the Chern number of a band j is 

Cj = J dk x dk y P(k x ,ky). (1.5) 

Here the integral is over the Brillouin zone. The total quantum 
Hall effect is then 


Oxy — 


h 



( 1 . 6 ) 


where in insulators each band is either fully occupied or com¬ 
pletely empty. 

The idea of the intrinsic anomalous Hall contribution is that, 
when a band is partially filled as in a metal, there remains 
something of the quantum Hall physics: there is a contribution 
to a xy of the same form as the Chern number, except that the 
integral is only over the occupied states: 

< rinSiC = J ^ J dk X dky F{k x , ky). (1.7) 

(For conciseness, here and below we simplify the equations by 
assuming one partially filled band and ignore the quantized con¬ 
tribution from other bands that are completely filled.) This 
integral can be transformed to a Fermi surface integral if de¬ 
sired 10 : We will review some of the standard arguments for 
this contribution in Section II. The complications arise because 
of the disorder required to achieve standard metallic transport 
(as opposed to Bloch oscillations, where the electron distribu¬ 
tion move periodically through the Brillouin zone with no relax¬ 
ation) . With disorder, there are additional contributions to a xy 
from impurity scattering (“skew-scattering” and “side-jump” 
are two mechanisms that scale differently with the diagonal 
conductivity. Disorder is particularly challenging to handle 
analytically in the two-dimensional case that we start with, as 
it is known that in a time-reversal breaking system (with or 
without spin-orbit coupling) the ultimate fate of noninteracting 
electrons at zero temperature and in the thermodynamic limit 
is always an Anderson insulator or a quantum Hall state. 

For this reason we would like to find a fully quantum treat¬ 
ment of anomalous Hall transport that is non-perturbative in 
the disorder strength and hence does not use the semiclassical 
equations of motion. It is known that if the disorder strength is 
kept constant as the system increases in size, eventually “quan¬ 
tum Hall plateau transitions” develop- 11 ^. (At very large dis¬ 
order, all states become localized 13 and there is only the ordi¬ 
nary ins ulat or.) The renormalization-group description of this 
procesJkEH leads to a two-parameter flow diagram, in which 
metallic conductivity appears only as an intermediate-coupling, 


nonuniversal effect. The Hall conductance of a finite-size dis¬ 
ordered system can be defined, under some assumptions, by 
an integral over periodic boundary condition^ 1 ®!, as reviewed in 
Section III. In the simplest case of a transition from ordinary 
insulator to n = 1 quantum Hall effect, virtually all realizations 
of disorder in a large system have a xy = 0 for filling less than 
some critical value n c and have a xy = 1 above this filling. An 
alternate case of Dirac fermions with spin-flip scattering where 
the Chern numbers move to band edges for strong scattering has 
been discussed^. A metal does appear in some other symmetry 
classes: for example, in the symplectic ensemble describing 2D 
systems with spin-orbit coupling, there is a metallic phase be¬ 
tween the ordinary and topological insulators (the latter is the 
quantum spin Hall phase®" 20 ). This phase shows up as a region 
of fillings over which the Ch ern n umber continues to fluctuate 
as the system becomes larger 

We will focus in this paper on approaching non-interacting 
metals in the unitary symmetry class (i.e., with broken time- 
reversal symmetry), essentially by scaling disorder to zero as 
system size grows. The details of this approach are provided in 
Section III. The basic idea is that, in addition to the plateau 
transition limit of strong disorder (more precisely, fixed disor¬ 
der as system size grows), there is a weak-disorder limit where 
Chern numbers are defined but fluctuate strongly from realiza¬ 
tion to realization. In other words, the Hall effect is strongly 
fluctuating between disorder realizations on a finite system. 
However, the average remains meaningful and shows nontriv¬ 
ial structure very similar to that of the metallic anomalous Hall 
effect. Section II reviews existing theories of the anomalous Hall 
effect and Section III introduces the theory of the Hall effect of 
a finite insulating system. There is a natural conjecture for the 
disorder-averaged Hall effect of insulators, which we compare 
to numerical simulations of a variety of models in Section IV. 
Section V discusses using the numerical results on the distribu¬ 
tion of pumping conductance as input to macroscopic models 
of conduction in disordered magnetic materials. An Appendix 
defines more precisely the mathematical question of how Chern 
numbers fluctuate between realizations and reviews some geo¬ 
metric considerations that enter into a simple model of the full 
probability distribution. 


II. KEY FEATURES OF THE ANOMALOUS HALL 
EFFECT 


The simplest approach to the anomalous Hall effect in met¬ 
als starts from the semiclassical equations of motion for a 
wavepacket of Bloch states. Remarkably, the AHE originates 
from a term in the semiclassical equations of motion that is ne¬ 
glected in almost all textbooks. This term was first obtained 
by Karplus and Luttinger, but partly because their work took 
place before the understanding of Berry phases, their results 
were not universally accepted. Consider the standard equations 

of motion 2 ^] 


Iik = eE + ev x B 

Hv = V fc e n (k) +- (II.l) 
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where ... indicate the often omitted Karplus-Luttinger term 
that we now explain. As an electron wavepacket moves in k- 
space under the influence of an applied field, there are two con¬ 
tributions to its spatial velocity v. The first term, the group 
velocity, describes how the modified energy-momentum relation 
changes the velocity of the center of a wavepacket, which is 
an effect that would be present even for a point particle. The 
Karplus-L utti nger contribution, which can be derived quite sys¬ 
tematically 1 I^1, describes how a change in k induces a change in 
the real-space location because the Bloch states are changing: 
this change is 

hv kl = — hit x .F(k). (II.2) 

Here T is the Berry curvature defined in the introduction, and 
for two-dimensional systems there is a single nonzero compo¬ 
nent of the Berry curvature. We will concentrate on the stan¬ 
dard idealized limit where the time-reversal breaking effect on 
observables is not through an orbital B held but rather through 
the Berry curvature JF. 

Now the simple, if not entirely convincing, approach to the 
AHE in the DC limit is to compute the current from the anoma¬ 
lous velocity term using just the applied held E without con¬ 
sidering the helds of impurities that normally act to produce 
a hnite conductivity. This could be justihed by applying the 
Kubo formula with an order of limits that is not quite what 
takes place in an actual DC experiment: if the applied held E 
is extremely weak, then it takes a long time for the electrons 
to accelerate significantly from their ground-state distribution, 
and during this initial window the transverse current will indeed 
be given by 


DC conductivity in some magnetic metals. For example, the 
orientation dependence of the conductivity tensor in anisotropic 
crystals is consistent with predictions from electronic structure 
calculations of the Berry curvatur^. 

One reason for taking the approach in this paper is to see how 
the AHE formula can be understood while avoiding the criti¬ 
cisms of other approached 27 . Part of this criticism in recent 
years grew out of consideration of the metallic spin Hall effectP^, 
where at least for some disorder potentials it appears that the 
Berry-phase contribution is cancelled to CP. The conclusion of 
the approach defined in the next section is that, for large finite 
samples, almost all of the variation of pumping conductance 
a xy with energy is explained by the intrinsic contribution. The 
only scope for a cancellation of the intrinsic part in the mea¬ 
sured conductance is thus if the Fermi-surface effects, which 
are expected to depend on details of scattering, exactly cancel 
the scattering-independent intrinsic part. The analytical pic¬ 
ture and numerics presented here thus support the viewpoint of 
Ref. [5] that, even in more realistic 3D models (note that again 
treating finite systems a 3D model becomes equivalent to a large 
number of bands in a 2D model) at nonzero temperature in the 
thermodynamic limit, the intrinsic contribution can be domi¬ 
nant. A way to state the resulting picture is that the pumping 
conductance computed below is more likely to be controlled by 
the intrinsic Hall effect contribution, as there is no Fermi surface 
so the other, scattering-specific contributions are minimized (or 
even zero)P 

III. PUMPING CONDUCTANCE AND STATISTICS 
OF CHERN NUMBER 


Jy = e 


d z k _ e 2 

Vy ~ 2t tH Ex [ 2tt 


(2tt) 


d 2 kT z 


(H.3) 


which is the same as Eq. 1.7 The integration is over the ground- 
state electron distribution, and time-reversal symmetry must be 
broken to give a nonzero answer since otherwise .F(k) is odd 
and E(k) even under k -A — k. An objection, possibly first 
raised by SrnitPS is that, in the usual picture of DC transport 
in a metal, the system has reached a non-equilbrium steady 
state (at least to linear order, i.e., neglecting Joule heating) 
between the applied electric fields and electric fields from static 
impurities. Then 7c is zero and the anomalous velocity term 
does not contribute. 

Hence the semiclassical approach might be useful at nonzero 
frequency but requires care at zero frequency, since a metal re¬ 
quires some level of disorder to have a finite DC conductivity. 
In other words, it is not a priori clear why it should be allowed 
to regulate the divergent diagonal conductivity by introducing 
a relaxation time r, while assuming that the relaxation’s effect 
on <j xy is controlled. Without disorder one would obtain Bloch 
oscillations, with no relaxation to equilibrium so that the usual 
assumption of the Kubo formula is invalid, and in the Bloch 
oscillation regime there would be no justification for just inte¬ 
grating over the ground state. Diagrammatic treatments can 
be quite complicated and do not seem to be in clear agreement. 
There is substantial theoretical and experimental evidenc^^ 


that the intrinsic formula (1.71 does explain a large part of the 


A. Hall conductance in mesoscopic systems via pumping 


In the thermodynamic limit of a two-dimensional 
independent-electron system with broken time-reversal 
symmetry, there is no metallic region of energies. Instead 
there are localized regions with occasional isolated energies 
supporting extended states of no nzero Chern number, and 
there is strong numerical evidence 111 3 that the localization 
length diverges near these critical energies as 
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1 

(e - E c y • 


(III.1) 


For a model system consisting of a tight-binding lattice Hamil¬ 
tonian with “topological bands” having nonzero Chern num¬ 
ber, to which on-site disorder is added, we expect at least two 
regimes. For very strong disorder, all states are localized and 
have Chern number zero, which could be viewed as happening 
via “floating” of extended states of opposite Chern numbers 
until they annihiliate 13 (note that the total Chern number of 
all the bands in a tight-binding model is always zero 30 ). For 
moderate disorder, as the system size grows larger with disor¬ 
der strength fixed, at least one sharp “plateau transition” en¬ 
ergy becomes well-defined within a topological band, and the 
collapse of the transition width (i.e., the range in energy over 
which Chern numbers flucutate) with system size can be used 
as a measurement of the localization length exponent v. 
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Both of these limits have only insulating phases and cease 
to show strong fluctuations of the Chern number between real¬ 
izations, except in a tiny range of energies around the plateau 
transitions. The focus of this paper is on the fluctuation regime, 
which can be obtained in a large system by running the disorder 
to zero as the system size increases. In order to explain the dom¬ 
inant effect observed in numerics in Section IV, we show how a 
simple theoretical argument leads to a conjecture for the aver¬ 
age Chern number in the fluctuating regime that corresponds 
to the intrinsic anomalous Hall formula discussed in Section II. 
In this section we give the theoretical motivation for the conjec¬ 
ture and discuss pumping versus Chern number. Fluctuations 
of topological properties between disorder realizations are well 
known to appear in numerical app roac hes to metallic regimes, 
and continue to be studied activeljEES b u t we are not aware of 
a previous quantitative analysis of the low-disorder limit. 

First, it is important to define what is meant by the zero- 
temperature conductance of a finite system. The quantity we 
co mpu te is the pumping conductance, as defined by Thouless et 
al.iaiLSJ f orma ]j ze ideas of Laughlin^! shown in Fig. [lj There 
are no leads and the transport is purely adiabatic: increasing a 
flux through the center of a Corbino disk, or through one great 
circle of a torus, leads to charge transport through the evolution 
of the ground state rather than through excitations. Adiabatic 
transport is possible because we expect no level crossings for 
a 2D system in the unitary ensemble because a degeneracy has 
codimension 3 (i.e., requires tuning three parameters) and there 
are only two free parameters (the boundary fluxes through the 
two great circles of the torus). 

The Landauer conductance for an open finite system with 
leads corresponds to a different physical picture: the levels of 
the finite system are broadened and there is a nonzero diago¬ 
nal conductivity, as opposed to the zero diagonal conductivity 
in the pumping picture. The way these two pictures become 
equivalent in the thermodynamic limit is that it ceases to be 
possible to remain adiabatic at any fixed nonzero frequency, as 
the energy scale of avoided crossings goes to zero. (Conversely, 
we expect the pumping conductance to be relevant to “closed” 
mesoscopic samples where the lead-induced linewidth is small 
compared to the level spacing.) The quantization of transport 
in the pumping limit means that we can analyze a xy in each 
disorder realization using the standard topological invariant, 
Chern number. The disorder-averaged pumping conductance 
then becomes a statistical average of different Chern numbers 
in different realizations. 


B. The intrinsic anomalous Hall formula as a statement 
about average pumping conductance 

The simplest way to imagine how the Chern numbers might 
be distributed under the influence of disorder is as follows. Re¬ 
turn to a lattice system with a topological band, say with Chern 
number +1, and for simplicity let us assume that this band is 
the lowest in energy. To study disorder, form a “supercell” of 
some large number N of the original unit cells, and now al¬ 
low on-site disorder (most kinds of local disorder that break all 
symmetries should behave similarly). Before adding disorder, 


“Pumping” conductance “Scattering” conductance 




Adiabatic (if non-degenerate) Dissipative (metallic) 

Closed system Open system 

No Fermi surface Fermi surface contributions 

States are evolving States are fixed 


FIG. 1: Comparison of pumping and scattering approaches to Hall 
conductance of a mesoscopic (i.e., finite) system. 


there would be N connected bands obtained from the original 
band, and the only well-defined topological invariant is the to¬ 
tal Chern number, still +1. Now by the argument above, there 
will be no degeneracies in the new Brillouin zone, which is N 
times smaller than the original unit cell but has N times as 
many bands. More precisely, the line crossings that are present 
in a supercell without disorder (since generically two surfaces 
will have linelike intersections) will no longer be present. In the 
absence of intersections, each of the N “minibands” will have 
its own Chern number, and all these Chern numbers will sum 
to Tl. 

At densities of k/N electrons in the lowest band that fill an 
integer number k of minibands in the Brillouin zone for the 
supercell, the pumping conductance will be given by the total 
Chern number of all filled minibands. So we are led to con¬ 
sider the statistical problem of how minibands pick up Chern 
number in different disorder realizations. The case of fixed dis¬ 
order strength in the thermodynamic limit is discussed above: 
for disorder not too strong, there is a well-defined plateau tran¬ 
sition at some density, and the metallic region of fluctuating 
Chern numbers around this transition decreases as a power-law 
in system size. 

When disorder is very small, so that its only effect is to avoid 
band intersections, the metallic region of fluctuating Chern 
numbers extends over the whole band, and both positive and 
negative miniband Chern numbers are occasionally seen. Con¬ 
sider two adjacent minibands. If on average the weak disor¬ 
der, in the process of eliminating band touchings, does not shift 
Chern number up or down in energy, then the mean Chern num¬ 
ber assigned to a band can come only from the Berry curvature 
for the corresponding density range in the non-disordered prob¬ 
lem. In other words, the mean Chern number for each mini¬ 
band will go to zero as ~ 1/iV, and the proportionality con¬ 
stant is related to the non-disordered Berry curvature for the 
corresponding interval of density (or presumably energy in the 
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weak-disorder limit). 

Hence, writing C/ for the Chern number of miniband i. this 
argument leads to the expression 


n F 


£<«> 


(a xy (E F )) = — Urn 

n F /N — n(E jp ) %— 1 

e 2 f F(k x ,k y ) 


'E{k x ,k y )<E F 


2?r 


dk x dky. (III.2) 


So we obtain the conjecture that the disorder-averaged quan¬ 
tized pumping conductance of an ensemble of weakly disordered 
insulators is far from being quantized and actually given by 
the (non-disordered) intrinsic anomalous Hall effect expression. 
The intrinsic anomalous Hall effect in metals is related to the 
averaged quantum Hall effect of insulators: this statement and 
its comparison to numerics are main results of the present work. 
Construction of a possible model for the full distribution of the 
pumping conductance at each filling, not just its mean value, is 
described in the Appendix. 

Clearly there are a number of assumptions in this argument 
that are difficult to justify: there is an assumption that the av¬ 
erage over weak disorder does not move Chern number in a pre¬ 
ferred direction in filling, for example. We will see in the follow¬ 
ing section, using recently developed numerical methods, that 
this prediction for the mean does correctly predict the dominant 
structure of the pumped conductance. An interesting question 
to which we do not have a full answer is about the full distribu¬ 
tion of miniband Chern numbers, beyond just the mean. Both 
positive and negative integers appear, and a first approach is 
to construct a model for the Chern number based on (a) the 
conjecture for the mean above, and (b) random resolutions of 
the crossings that make the miniband Chern numbers indepen¬ 
dent and are almost statistically independent of each other (see 
Appendix). Since the simplest experimental observable by far 
is the mean conductance, we will focus on numerical tests of the 


conjecture (III.21 here. 



FIG. 2: Intrinsic contribution to the Hall conductivity in the 1/5 
Hofstadter problem as a function of filling fraction, calculated from 
an integral over the Berry curvature, Eq. (IV.9I. Dotted lines indi¬ 
cate when a band is completely filled and a new band begins. Inset: 
band structure of the model. 


To study larger system sizes and obtain a reasonable test of 
the conjecture regarding Chern numbers, we use the method of 
Hastings and Lorin^lM) Let p again be the projection operator 
onto some set of bands. Then define two matrices U and V using 
the band-projected position operators in the two directions: 

Pexp(*0)P ~ ^ , Pexp(*$)P~ (IV.2) 


An integer (believed in general to equal the Chern number) 
can be defined that characterizes how far the U and V matri¬ 
ces, which are almost unitary and almost commuting, are from 
commuting: 


Tr log (VUV^U^) = 27 rim + r, meZ (IV.3) 


IV. NUMERICAL METHODS AND RESULTS 


or alternately 


A. Approaches to efficient computation of Chern number 
for pumping conductance 


det(VUV^U^) = exp(27r im + r), to € Z 


(IV.4) 


The first approach used to calculate Chern number was the 
conventional one of an integral over momenta (equivalently, 
over boundary phases® associated with phase-periodic bound¬ 
ary conditions on the torus). In the supercell, we define the 
projection operator P = \u l k )(u l k \ onto the subspace of miniband 
i. Following RefF®. we can then compute the Chern number as 
the integral over momenta 

C z =- l -J dk x dky Tr P(d x Pd y P - d y Pd x P). (IV.l) 

With a sufficiently fine mesh, this approach indeed gives quan¬ 
tization of subbands, but the fine mesh required means a con¬ 
siderable computational effort, limiting accessible system sizes 
to 10 by 10 or less. 


where the definition of the determinant is by exponentiaing the 
trace of the logarithm as in the previous equation. The integer 
to can be nonzero once we note that there are many eigenvalues 
so that the total of their imaginary parts may be nonzero and 
equal to 2-Km even if each individual eigenvalue is . 

A feature of this approach is that the method is constrained 
to return an integer for each band, even in the clean case where 
no Chern numbers for individual minibands are defined. Hence 
we expect that at sufficiently small disorder the results cannot 
be trusted, and indeed at disorder values smaller than the ones 
in the following section, the Chern numbers become seemingly 
random small integers. We have restricted ourself in the follow¬ 
ing to large enough disorder that the Chern number results are 
not distorted by this type of systematic error. 
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B. Tight-binding models with nonzero Berry curvature 


The first model we use to illustrate our findings is the typical 
Hofstadter problem on a square lattice, 


H i = 


H [ c j+i k c ik + e ijAo c] k+1 c jk + h.c. 


(IV.5) 


j,k 


in terms of creation operators for particles on site (j-k) 
of the square lattice. Electrons can hop from site to site while 
motion in y-direction comes hand in hand with a site-dependent 
phase e _IJj4 ° that reflects the presence of a uniform magnetic 
field perpendicular to the lattice plane. We concentrate on the 
special value Aq = ( f > o/5 of the magnetic field, where each unit 
square contains one-fifth of a flux quantum and a unit cell of 
five sites can be defined. The electronic behaviour can then be 
described by a band structure with five bands as shown in the 
inset of Fig. [2j The figure itself shows the intrinsic contribution 
to the Hall conductivity as a function of the filling fraction of 
the bands, following the Berry curvature integral in Eq. (IV.91. 


For filled bands (dotted lines), it equals the sum over the Chern 
numbers of the occupied bands. In this model, the bands carry 
Chern numbers 1,1, —4,1 and 1 respectively. 

The second model we consider is again a hopping model on a 
square lattice but this time with a spin degree of freedom. It is 
best defined for a clean system in momentum basis first, 


H 2 (k) 


n(k) ■ <x = 


( sin k x 
sin k y 

m — cos k x — cos k 



with Pauli matrices cq and momenta k x ,k y £ [— n, n) in the first 
Brillouin zone of the two-dimensional system. The two bands 
of the model, E±(k) = ±|n(fe)|, are plotted on the left in Fig. 
[3] for m = 1.9. 

The corresponding hopping model reads 


^2 — 2 [ m £scjfc s Cjfcs ^,s[Cj +lks Cjks + Cjk+ls C i ks \ 


j,k,s 


* C j+l ks C jks + £sCj k + 1 gCjks + h.c. (IV. 7) 


in terms of creation operators c] fcs for particles on site (j,k) of 
the square lattice with spin s. Here = 1 and £ 4 . = — 1 help to 
distinguish terms with different spin. 

While the hopping model may seem complicated, the specific 
beauty and usefulness of this model comes with the fact that 
only two bands exist and that the Berry curvature can be ex¬ 
pressed in terms of the normalized vector n(k ) = n(fe)/|n(fe)|, 
yielding 


a 


intrinsic 

xy 



occ 


dh(k) 

dk x 


x 


dn(k) 

dk y 


■n(k), (IV. 8 ) 


when only states in the lower band are occupied. The right 
part of Fig. [3] shows the intrinsic contribution as a function of 
filling for this clean hopping model. For a completely filled lower 


band, the intrinsic contribution cr^ / trinslc = le 2 /h is determined 
by the Chern number 

I(ra ) = {s, g; ( m) ifH<2; (IV . 9) 

The vector h(k) can be interpreted as the spin quantization 
axis for a given momentum and the Chern number becomes the 
winding number of this vector throughout the Brillouin zone, 
as illustrated in the inset of Fig. [3] for m = 1, where 1=1. 


C. Numerical results with disorder 

Disorder is introduced in our models by adding a random 
potential energy on each site of the lattice to the Hamiltonian, 


— H\ + Y Ej k c\ k Cj k , 
j,k 

(IV.10) 

H 2 = H 2 + Ej ks c*j ks Cjks- 

(IV.ll) 


(IV.12) 


We take these random energies from a gaussian distribution 
with zero mean and variance U. One can reduce undesired lo¬ 
calization effects slightly by considering a different type of dis¬ 
order distribution. We have tested this for a Wigner semicircle 
distribution - the results did not change much though. Our 
presentation will thus focus on Gaussian disorder. 

To confirm our hypothesis numerically, we take disordered 
square samples of size L x L for either model and calculate the 
Chern numbers for the disorder bands following the methods 
mentioned in Sec. |IV A[ We then plot the disorder average of the 
sum of Chern numbers of the first l disorder bands against the 
“filling” n = l/L 2 and compare it to the intrinsic contribution 
to the Hall conductance of the clean case. 

Generally the numerical treatment has to be done carefully: a 
certain size is required so that sufficiently many disorder bands 
guarantee applicability of a statistical description. On the other 
hand Anderson localization sets in at a certain sample size for 
fixed disorder strength. This has to be avoided to maintain a 
metallic behaviour. Furthermore the numerical procedures only 
work reliably for large enough disorder strengths. For reasons of 
numerical efficiency, all results shown below have been obtained 
with the method of Hastings and Lorin^H _ agreement with the 
conventional method of integration over boundary fluxes has 
been established for smaller sample sizes though. 

Fig. a shows our results for the disordered one fifth Hofs¬ 
tadter problem. We have focussed on the first two bands where 
all general trends can be observed. For the first band, the we 
find remarkable agreement between the disordered data and the 
intrinsic Hall conductivity in the clean case (red line) at sam¬ 
ple sizes L = 20,30,40 and U = 0.01. One can see that the 
agreement becomes better with growing sample size. For the 
second band, the deviations between the clean and disordered 
cases are still a bit larger although the shape shows a clear 
resemblance. Naturally different bands may show a varying de¬ 
gree of convergence between the two cases. We suspect that 
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FIG. 3: Left: band structure of the two band model in Eq. (IV .61 for m = 1.9. Right: intrinsic contribution to the Hall conductivity as a 
function of filling fraction of the lower band as calculated from an integral over the Berry curvature, Eq. (IV.91. As one can see at integer 
filling, n = 1.0, the band carries a Chern number of one for m < 2.0, while above the topological phase transition at m = 2.0, the Chern 
number becomes zero. Inset: Spin quantization axis of the model in momentum space, illustrating the non-trivial winding that exists when 
the Chern number is one (m = 1.0). 


larger samples (paired with smaller disorder) would bring the 
disordered curves closer to the red line - this is presently be¬ 
yond the power of our numerical treatment. One can further 
see some bumps in the green (L = 20 ) curve in the second band 
- these are artifacts of too small sample sizes where the sta¬ 
tistical treatment breaks down (the disorder band bandwidths 
are too large and the bands overlap strongly). For even smaller 
samples (e.g. L = 10), the curves deviate strongly and show 
additional undesired signatures. 


The lower plot of Fig. [4]compares different disorder strengths. 
Disorder strengths U = 0.01 and U = 0.001 give practically the 
same results, the latter may already be effected slightly by nu¬ 
merical imprecisions. A deviation can be seen for U = 0.1. In 
the lower band, localization sets in and flattens out the contri¬ 
bution of the outer disorder bands - around Filing 0 and 1. The 
accumulation of Berry curvature is then concentrated on the 
band center. When investigating the behaviour of the second 
band, it may surprise to see better results for larger disorder. 
While the opposite behaviour would have been more reassur¬ 
ing, this doesn’t constitute a problem per se, since there is an 
easy explanation: our conjecture is expected to hold for the 
small disorder - large system limit where the length of the sys¬ 
tem needs to be rescaled by L oc 1/U 2 with decreasing disorder 
strength for curves to be comparable. Thus curves at fixed dis¬ 
order are not required to show any clear relative behaviour when 
compared. An optimal treatment would show curves for differ¬ 
ent disorder strengths and fixed U 2 L = const. Unfortunately 
the rather limited numerically accessible range of system sizes 
(L < 40) renders such an approach fruitless. 


Let us now proceed to the second model (Eq. IV. 6 |. Fig. [5] 
demonstrates the good agreement we obtain between the dis¬ 
ordered data and the intrinsic Hall conductivity in the clean 
case (red line) at various sample sizes and disorder strentghs. 
We chose m = 1.0 for the plots since then the gap between the 
bands is relatively large and the bands carry Chern numbers 
1 and —1 respectively. Even though it is hard to see, larger 
samples provide better agreement (upper plot), just as lower 


disorder strength does. The green curve in the lower plot shows 
deviations due to localization of the outermost bands (around 
Filing one), caused by the strong disorder U = 0.1. The two 
band model further allows us to explore the situation where the 
chern number of the bands is zero but where there are still Berry 
curvature effects. In Fig. [ 6 ] we show data for m = 2.2 where 
this is the case. Once again we Fnd reasonable agreement which 
grows with system size. 

In summary our numerical results support our theoretical pre¬ 
dictions qualitatively. A full quantitative agreement is out of 
reach due to the computational costs required to reach full con¬ 
vergence. 


V. METHODS TO MODEL EXPERIMENTS 

In this section we discuss briefly how the purely quantum, 
T = 0, Fnite-size results of this paper could be used as a basis 
to model experimental parameters under more realistic condi¬ 
tions. At T = 0, with fully coherent transport, there is no metal 
in two dimensions with broken time-reversal symmetry; rather, 
in the weak disorder limit discussed here, there would be strong 
sample-to-sample fluctuations between ordinary and quantum 
Hall insulators. At nonzero temperature decoherence leads to 
a nonzero conductivity. The problem of obtaining the effective 
conductivity tensor of disordered materials under time-reversal 
symmetry breaking is a widely studied one. These studies typ¬ 
ically fall into two categories - one that studies a system with a 
continuously varying conductivity tensor, and one that refers to 
a system that is broken into regions with different conductivi¬ 
ties and a sharp boundary between them, a so-called composite 
medium. The result of the preceding sections, that a Fnite-size 
region of a disordered sample has a statistical distribution of 
pumping conductance in the weak-disorder limit, can serve as 
inputs to these larger-scale models. 

The effective conductivity tensor of a system is deFned as 

O') = °e(E) 













0.0 L^T.. 1 

0.0 0.5 1.0 1.5 2.0 


n 



FIG. 4: Disorder-averaged Chern sum of the first l disorder bands 
plotted against ’’filling” n = l/L 2 for square samples of size L x L 
in the one fifth Hofstadter problem. Upper: sample size dependence 
for disorder strength U = 0.01. Lower: disorder strength dependence 
for L = 40. Averages were taken over more than 2500 different com¬ 
binations of random disorder configurations and random boundary 
phases. The red lines represent the intrinsic contribution to the Hall 
conductance in the clean case. 


where (j) and ( E ) are the spatial (or volume ) averages of the 
local current density, j(r), and electric field, E(r). These in 
turn are connected through the local relation j(r) = <r(r)E(r ), 
with cr(r) the local conductivity tensor. The first approach 
to obtainin g cr e is known as the effective medium approxima¬ 
tion (EMApll25J similar in spirit to a mean-field approxima¬ 
tion, where a composite system is thought of as being made 
of grains with local conductivity cq, embedded in a uniform 
background material with conductivity a e , which is solved for 
self-consistently. This method is most useful when considering 
a system made of a random mixture of two phases with scalar 
conductivities, where a closed form solution can be obtained for 
a e as a function of the two local conductivities and the fraction 
of space they occupy. For non-scalar conductivities, which is 
the case of interest for systems with broken time reversal sym¬ 
metry and a Hall-like response, this approximation allows for 
asymptotic solution at very high values of the magnetic fielcP®. 

Other methods focus on studying inhomogeneous quantum 
Hall systems near a plateau-plateau transition, taking a perco¬ 
lation theory approach where the picture of current networks 
forming at the boundary between quantum Hall puddles of dif¬ 
ferent filling fraction is used to predict the scaling behavior of 
the Hall resistance!- 37 ^ Two dimensional inhomogeneous media 



n 



n 


FIG. 5: Disorder-averaged Chern sum of the first l disorder bands 
plotted against ’’filling” n = l/L 2 for square samples of size L x 
L in the spinful two band problem at m = 1.0. Upper: sample 
size dependence for disorder strength U = 0.001. Lower: disorder 
strength dependence for L = 20. Averages were taken over more than 
10000 different combinations of random disorder configurations and 
boundary phases. The red lines represent the intrinsic contribution 
to the Hall conductance in the clean case. 


in a magnetic field can also be studied through using linear frac¬ 
tional transformations and duality relations 35 . For example, in 
quantum Hall systems, fluctuations in the background potential 
cause the system to phase-separate into puddles of two quan¬ 
tum Hall phases with different Hall conductance. One way of 
applying linear fractional transformations for this scenario is by 
thinking of the local conductance as a complex number rather 
than a tensor, eq = o’ , xx +ia l xy and mapping the system onto one 
with with zero magnetic field, where for a two phase compos¬ 
ite exact solutions for the effective conductivity can be found. 
This way one can also derive a semi-circle relation between the 
longitudinal and the Hall conductance (another way is to map 
the system onto itself but flipping the direction of the magnetic 
fieidpUMHI] 

For more complicated scenarios where we expect the system 
to be composed of a random mixture of phases obtaining multi¬ 
ple values of the conductivity, a network model was formulated 
and numerically solved by Parish and LittlewoocP^ originally as 
means to account for the non -saturating linear magnetoresis¬ 
tance observed in some inhomogeneous semiconductors. Within 
their model, the material is discretized into a random resistor 
network in a perpendicular magnetic field, where each unit is 
a four terminal disk. Four incoming and outgoing currents, Ii 
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FIG. 6 : Disorder-averaged Chern sum of the first l disorder bands 
plotted against ’’filling” n = l/L 2 for square samples of size L x 
L in the spinful two band problem at m = 2.2. Upper: sample 
size dependence for disorder strength U = 0.001. Lower: disorder 
strength dependence for L = 20. Averages were taken over more than 
10000 different combinations of random disorder configurations and 
boundary phases. The red line represents the intrinsic contribution 
to the Hall conductance in the clean case. 

are connected to four voltage differences Vj via an impedance 
matrix Z, 

Vi = Zijlji 

which is found by solving the Laplace equation with the appro¬ 
priate boundary conditions. The parameters in Z vary from 
unit to unit, generating the randomness of the network. The 
overall impedance of the network is the found numerically by 
using Kirchoff’s law. The interesting behavior observed when 
mapping out the local current profile is that current loops form 
in the bulk, producing a Hall-like behavior with a Hall resistance 
linear in the magnetic field in the infinite network limit. 


VI. CONCLUSIONS AND FUTURE DIRECTIONS 

This paper demonstrates a close connection between disorder- 
averaged pumping conductance in weakly disordered mesoscopic 
systems and the instrinsic anomalous Hall effect in metals. Only 
greatly simplified model band structures have been studied, and 
an obvious direction for future work is to see if the same strong 
correlation between pumping conductance and the intrinsic Hall 
effect formula is observed in realistic models of magnetic ma¬ 
terials such as iron. On the theoretical side, the Chern num¬ 
ber is just one case of a topological invariant that is additive 
over bands, and the same question of how the topology is ap¬ 
portioned over minibands can arise for other symmetry classes 
and topological invariants. A specific example is the symplectic 
case (i.e., a material with time-reversal symmetry and spin-orbit 
coupling), which in two dimensions has a Z 2 invariant. If the 
spin-orbit coupling conserves a U( 1) symmetry such as spin in 
the z-direction, then the problem will reduce to two copies of 
that studied in this paper. For generic spin-orbit coupling, the 
statistical properties of Z 2 invariants are an interesting question 
that requires a different theoretical approach than applied here. 


The improvement of numerical methods 33 ’ to calculate topo¬ 
logical invariants on a large system was an important enabling 
development for the present work, and we hope the present ap¬ 
plication will stimulate further methodological developments. 
Another important theoretical development was the recent suc¬ 
cessful analysis of renormalization-group flow of conductivity 
in universality classes AIII and BDI in one dimension^®, which 
have topology similar to that of the integer quantum Hall effect. 
A similar analysis for the quantum Hall universality class should 
yield a limit of non-quantized a xy obtained by reducing disor¬ 
der with increasing system size as done here. It would be quite 
desirable to develop a field-theoretic treatment of disorder aver¬ 
aging in the metallic regime that incorporates the microscopic 
information contained in the Berry curvature on momentum 
space in order to reproduce the observed dependence of a xy on 
density. 

The authors acknowledge useful conversations with A. Altand 
and T. Loring and support from a Rubicon fellowship (J.D.), 
AFOSR MURI (R.I.) and grant NSF DMR-1206515 (J.E.M.). 
Parts of this work were completed at the Moorea Center for 
Advanced Studies. 


Appendix A: Geometrical considerations and a model for 
statistical distribution of Chern number 

Note : This mathematical appendix was written by D. Freed 
and M. Freedman. It provides a more precise mathematical 
definition of the main question defined in the main text and 
gives some geometrical background that is useful in understand¬ 
ing how Chern numbers are statistically distributed. A simple 
model to capture these considerations is then defined. 

As in the main text consider a two-dimensional planar sys¬ 
tem H symmetric under a discrete lattice Z ® Z of translations 
with Brillouin zone (= momentum torus) T = Tbi g . Suppose, 
whimsically, we consider a larger “super” unit cell, SC , let us 
say consisting of n 2 original unit cells C. Correspondingly there 
is a small Brillouin zone T sm aii and a covering map 


Ismail 

where the n 2 -sheets of the cover are permuted by the action of 
the covering group G = fine lattice/course lattice = Z n © Z n . 
Each character (homomorphism to t/(l)) on the coarse lattice 
lifts to n 2 different characters on the fine lattice. 

We consider a single low energy band of the system Hamilto¬ 
nian H giving a state (eigenvector of H^) and an eigenvalue 
efe, k £ Tbig- Thus e is a function, e : Tbi g —» R, and determines 
an n 2 -sheeted “multi-function” e := e o 7 r _1 : T sma n —► R. 

The reason for considering the band as a multi-function on 
T sm aii is model disorder as: 

1. a passage from C to super cell SC, to reflect reduced 
symmetry, and 

2 . a resolution or “avoidance” of the crossings of e. 
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This resolution produces precisely n 2 -bands e*, i < 1 < n 2 , 
e-i '■ ?smaii —> R- For topological reasons (the additivity of Chern 
classes -®) the Chern number is conserved: 

n 2 

ci{e) ='^2c 1 {e i ). (A.l) 

i —1 


The ordered list of Chern numbers: (ci(ei), Ci(e 2 ),..., ci(e n 2 )), 
the “Chern sequence” c, while deterministic in the precise choice 
of disorder, can also be regarded as a random variable on a mea¬ 
sure space of possible disorder realizations (taken from disorder 
model in which parameters, e.g., “strength of disorder”, have 
been fixed). 

It looks exceedingly hard to solve a general instance of the 
deterministic problem so it may be of some value to propose 
a simple statistical model which captures key aspects of band 
resolution and requires as input only the most basic features 
of the band e and the disorder model. We write a Gibbs-type 
distribution for the Chern sequences c by assigning to each c 
a kind of mock energy E(c), with p(c) = ^e~^ E<F> , where the 
partition function Z = ')2 c e~^ E( ' c \ 

Strangely, high /? (i.e., low temperature) corresponds to 
strong disorder, because as explained in the main text large 
disorder leads to the standard quantum Hall plateau transition 
at a fixed density and energy. The fluctuating regime corre¬ 
sponding to high temperature is physically the metallic regime 
of weak disorder. 

As input to the model we need two functions on T 2 ig : the 
energy e : T 2 ig —> R introduced above and the Chern density 
c : Ibig —► R. The state ipk is a mapping ipk : T 2 ig —>• BU(1) into 
the classifying space consisting of rank = 1 Hermitian projectors 
ipk l— t Pk ’■= | ipk) (V’fe I- Concretely the Chern form, w, pulls back 
to '-L big as. 


w := -— tr ( P 
2m 


dP dP 
dk T dk„ 


dP dP 


dk y dk x 


dk x A dky , (A.2) 


which specializes to the integrand of line (20) of this paper when 
the projectors act within C 2 . Formula A.2 also gives the C\- 
integrand on all BU{m) for general rank. Besides e and c there 
is a parameter r £ R + U 0 that (roughly) measures influence of 
generic e-triple points: r increases if there are more k £ T^ ig 
with e{k) = e(g 1 {k)) = e{g 2 (k)) for k £ T£ g and g 1: g 2 ± id are 
distinct elements of G := Z n x Z n . The final input is a choice of 
reciprocal temperature /3 « disorder strength. (A more precise 
understanding of this relationship would be desirable.) 

Before describing the model we step back and give some ab¬ 
stract background for the phenomena we propose to model. The 
single band {ipk} defines a complex line bundle E —► T 2 ig , i.e., 
the line Ek = C ipk, the span of ipk , Fes over the point k. E 
pushes forward to a rank n 2 vector bundle F —► T 2 mall where 
the fiber F :] over j £ T 2 is: 


F :i ~ @ E ki 

he- 7r _1 (j) 


the direct sum being unordered. Each summand is still labeled 
by its energy eigenvalue e(k). 


Topologically F is isomorphic to a direct sum of n 2 line bun¬ 
dles. We give two reasons for this: the first abstract, the second 
concrete. 

First, the classifying space BU(n 2 ) in dimensions < 2 n 2 is 
homotopy equivalent to stable BU, which by Bott periodicity 
has: 


n 2(i-i){BU) = 0, n 2 i (BU) = Z, i> 1. 


(A.3) 


Thus any map Tbi g —> BU[n 2 } : for dimensional reasons, factors 
through the 2-skeleton of BU (n 2 ) which it shares with BU( 1): 

SI< 2 (BU(n 2 )) = SK 2 {BU( 1)) = CP 1 “ S 2 . 

This tells us F actually has the form (line bundle) © 
(trivial line bundles). This is more than we want. We want 
to arrive at a physically natural decomposition which in gen¬ 
eral will have many nontrivial factors. 

The concrete picture is based on “eigenvalue repulsion” which 
states that in families of Hermitian matrices collision of eigenval¬ 
ues is a codimension 3 event. Since dimension(T 2 mall ) = 2 < 3, 
the multi-bands of F, e, will not collide if allowed, through dis¬ 
order, to become generic. Then the eigenvectors associated to 
n 2 disjoint bands over T sma u define a decomposition into line 
bundles F = © line bundles. 

Let us look first at the topological picture of resolving a multi¬ 
band in one dimension, for, say, n = 3 (see Figure [7|. In two 
dimensions the situation is similar but there are now generically 
triple points of bands to consider. Looking down from the top, 
say, at the graph of the multi-function e : T 2 mall -)Kwe see 
triple points as in Figure [ 8 ] 

The resolution of triple points together with the previously 
shown one dimension resolution (Figure [ 7 ]), but now in 1- 
parameter families along double point arcs and double point 
circles, occurs when the Hamiltonian Hk is perturbed to avoid 
eigenvalue collision. The result is some explicit way of decom¬ 
posing F as a sum of line bundles. By (A.l) the total Chern 


number is conserved, and this is all we know rigorously. But 
we now present a detailed heuristic for how to understand this 
conservation law as locally implemented by the resolution of 
double and triple points. This perspective, although heuristic, 
suggests a way to model the “diffusion” of curvature between 
nearby small bands. 

In mathematics, symmetry often destroys genericity and so it 
is here: the G = Z n x Z n symmetry of the cover T 2 ig -£ T 2 maU 
forces a non-generic eigenvalue distribution on F —> T 2 mall . Dis¬ 
order breaks this symmetry and perturbs the n 2 bands to be 
disjoint, much as drawn in Figure [ 8 ] The weaker the disor¬ 
der the closer the classifying map comes to the co-dimension 
= 3 eigenvalue-collision locus X , so the more violently the 
eigenspaces must move in following the resolution which car¬ 
ries them between orthogonal eigendirections. If the resolution 
is on a (momentum space) length scale l. the velocities of the 
eigenvectors scale as P 1 . This leads to our interpretation of 
effective temperature in the model. 

A final preliminary to defining the model is to understand a 
conservation law pertinent to the resolution of a single double 
circle between two sheets of e. Consider an angular region A C 
^smaii around the double curve, and for simplicity assume no 
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FIG. 7: One-dimensional illustration of the map from big to small Brillouin zone; these correspond to single unit cell and supercell in real 
space. 



FIG. 8 : Resolution of triple points into disjoint disks under weak disorder. 


other double curves meet A. Imagine the disorder perturbs A 
relative its boundary dA. Before the perturbation interior A 
meets X (actually along a circle but this is not relevant). The 
universal first Chern form w for the two bands integrates to 
±1 respectively on the linking 2-sphere to X; it is the Bott 
generator. 

The perturbation to A relative DA will generically contribute 
oppositely: ± /y w to the two bands. The contribution is not 
generally integral since dA = dA 0. Thus when bands are re¬ 
solved we expect bits of curvature to be added to one band and 
subtracted from an adjacent band. Triple points make the pic¬ 
ture still more complicated by sharing curvature between triples 
of adjacent bands. We model this effect by introducing a corre¬ 
lation coefficient, r. We now define the model. 

Over T s 2 mall we define the pre-resolution bands e', 1 < i < n 2 , 
as e'(j) = i th largest value among e(k) where n(k) = j. The e' 
are piecewise smooth and we define a,; = f, w. Note that {aj | 

1 < i < n 2 } can be computed from the functions e and c. Now 
consider all possible “double point” corrections di, 1 < * < n 2 , 
so that the following numbers are integers: 

c = ci = ai + di, c 2 = a 2 - di + d 2 , c 3 = 03 — d 2 + d 3 ,..., 
Cn 2 = 0-n 2 - d n 2- 1 - (A.4) 

Clearly d\ € Z — a, d 2 € Z — a 2 + d\, ... The c’s represent the 
Chern classes on the resolved bands ej. 

The “mock energy” we use is: 

n 2 — l 

E{c) = X! + T \ d i ~ d i+ ll) 1 


where by convention d n 2 = 0 . 

Remarks. The philosophy is that the resolution as a map to the 
classifying space 

2 2 2 
n n n 

e = II e - II ->• l[BU(l) C BU(n 2 ) 

i—1 i— 1 i—1 

lias a geometrically natural energy, J d area | Ve| 2 , of which the 
portion integrated away from the multiple points is fixed by the 
data e and c. For the integral near the double points we use the 
Chern curvature density as a surrogate for energy. There is a 
Wirtinger-type inequality: 

|Chern curvature e| < const. |Ve | 2 , (A.5) 

and we are using the left-hand side to estimate the right-hand 
side. This is likely a better approximation the higher the dis¬ 
order. Disorder sets the length scale (in momentum space) 
over which e can relax toward a map using only the topologi¬ 
cally mandated energy. Similarly, the longer this length scale, 
the more the perturbed system explores its phase space before 
choosing the curvatures d. This is the reason increasing disorder 
corresponds to increasing /3. 

The t \d- t — dj i+ \ term is an attempt to capture the influ¬ 
ence of triple points, which spreads out curvature over adjacent 
triples of bands. The term enforces some local correlation be¬ 
tween the di and di + \. 

The appearance of absolute values may look strange but this 
is the correct scaling for energy/Chern number. To see this, 
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consider a composition 

f 1 * °') 

T 2 —i-T 2 -U- BC/(!) 



g*C\ [T 2 ] = n 2 f*Ci [T 2 ] and / T2 darea||Vg || 2 = 

n 2 J T 2 d area11V/11 as well. 

There is an invisible factor n 2 ~ 2 in the definition of E. Pass¬ 
ing from a band on Tj^ ig to a similar band on the smaller T 2 mall , 
the integrand of energy scales as n 2 but the integration measure 
as n~ 2 . 

Given E , define the clearly convergent partition function: 


Z(fi) = Y^c e ~^ E ^> where the constraint on c is ^ c = ci(e). 
Then the model probability for observing a vector c of Chern 
numbers for the n 2 bonds over T 2 mall is: 


p(c1 = 

As explained above, /3 ~ disorder strength, as at large dis¬ 
order the distribution ceases to fluctuate corresponding to low 
temperature. For a wildly fluctuating band b, r should be larger 
since there will be more triple points of e, but for ordinary band 
geometries, r may be fairly universal. 
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